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Abstract. We present a new method to construct an action functional for a field theory 
described in terms of nonlinear partial differential equations (PDEs) . The key idea relies on an 
intrinsic representation of the PDEs governing the physical system relatively to a diffeomorphic 
flow of coordinates which is assumed to be a functional of their solution. This flow, which will 
be called the conjugate flow of the theory, evolves in space and time similarly to a physical fluid 
5— I , flow of classical mechanics and it can be selected in order to symmetrize the Gateaux derivative 

£-4 of the field equations relatively to a suitable (advective) bilinear form. This is equivalent to 

require that the equations of motion of the field theory can be derived from a principle of 
stationary action on a Lie group manifold. By using a general operator framework, we obtain 
the determining equations of such symmetrizing manifold for a second-order nonlinear scalar 
field theory. The generalization to vectorial and tensorial theories is straightforward. 
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1. Introduction 

The problem of finding an action functional whose Euler-Lagrange equations correspond to a 
prescribed set of partial differential equations (PDEs) is known as the inverse problem of the 
calculus of variations and it has attracted the attention of researches for more than one century. 
Perhaps, one the main reasons is that the formulation of a field theory in terms of an action 
functional is very elegant and, more importantly, it allows us to establish an immediate connection 
between symmetry principles and conservation laws |231 132) . As is well known, the existence 
conditions of an action functional can be put in a correspondence with the theory of irrotational 
vector fields [5^1 SSI ■ Essentially, if the path integral of the nonlinear operator representing 
the field equations is independent on the trajectory of functions connecting two specified points 
in a function space (i.e. in the domain of the nonlinear operator), then there exist a scalar field, 
the action, whose variational derivative is stationary in correspondence of the field equations 
of the theory. The path integral of an operator along a trajectory of functions is defined in 
terms of a bilinear form [25j[4TJ. Thus, the solution to the inverse problem of the calculus of 
variations, i.e. how to determine the action functional of a given set of PDEs, can be reduced to 
looking for a bilinear form that makes a given nonlinear operator "irrotational" , if any. It has 
been shown by Tonti [42] that there exist not just one but an infinite number of bilinear forms, 
often depending on the field equations of the theory, that satisfy this requirement. Therefore, 
an infinite number of action functionals can be constructed for a given set of PDEs. However, 
the physical meaning of the corresponding action principle sometimes may be obscured by the 
generalized bilinear forms that have to be adjusted on the specific problem. An alternative way 
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to proceed is to select a specific bilinear form having a physical meaning, and then look for ways 
of modifying the given field equations as to obtain a new problem which is potential with respect 
to the chosen bilinear form. Among known methods devised to do so, we recall the method of 
adding the adjoint equations [30l [13 [23 EH] and the integrating operator method of Tonti [42] . 

The purpose of this paper is to introduce a new approach to construct an action functional 
for an arbitrary field described in terms of nonlinear partial differential equations . The key idea 
relies on an intrinsic representation of the field equations relatively to a flow which is assumed 
to be a functional of their solution. This flow will be called the conjugate flow of the theory. 
Let us briefly describe the main ideas that led us to introduce this conjugate flow and, more 
importantly, their relevance in the context of known physical systems. To this end, we first 
notice that flows of coordinates being functionals of the solution to a field equation or a system 
of field equations arise naturally in many areas of mathematical physics. Perhaps, the most 
relevant example is in the context of classical fluid mechanics, where the set of fluid element 
trajectories in space is related to the velocity field that solves, e.g., the Navier-Stokes equations 
[TJ [2j . This flow of curvilinear coordinates is known as Lagrangian coordinate system and it is 
obtained by integrating the definition of the velocity field [33]. In this sense, physical fluid flow 
of classical mechanics can be considered as a very particular type of conjugate flow. Another 
example is the free-falling coordinate system 49 in the Einstein's theory of gravitation. The 
conjugate flow here appears as a geodesic flow |23j in a four-dimensional Ricmannian manifold 
whose metric is determined by a particular distribution of energy and momentum through the 
Einstein's field equations. 

In both examples above, the relation between the solution to the field equations of the 
physical system and the flow of curvilinear coordinates is known and it reduces to the definition 
of the velocity field in the case of Navier-Stokes equations and to the definition of geodesies in 
the case of Einstein's theory of gravitation. In a broader framework, however, such functional 
relation may be left unspecified. This key observation provides us with an infinite number of 
functional degrees of freedom (those associated with the conjugate flow) that can be selected, 
e.g., by requiring that the field theory, expressed in conjugate flow intrinsic coordinates, is 
derivable from a principle of stationary action. In other words, we are posing the following 
fundamental question: does it exist a junctional flow oj curvilinear coordinates such that ij we 
represent the given set oj field equations relatively to that flow then an action functional can be 
constructed? This is equivalent to look for a representation of the field equations on a Riemannian 
manifold [7] 134) that depends on their solution, in such a way that the Gateaux derivative of 
the nonlinear operator associated with the field equations (in intrinsic coordinates) is symmetric. 
This formulation of the inverse problem of the calculus of variations brings together concepts 
of differential geometry and nonlinear functional analysis and, as we will see, it results in the 
formulation of new types oj action principles generalizing those ones based on specific functional 
flows of coordinates, such as the Herivel-Lin principle (T9l [TTl IB] . 

This paper is organized as follows. In section [5] we introduce the theory of the conjugate 
flow and we characterize the group of infinitesimal perturbations by using methods of nonlinear 
functional analysis. The representation of arbitrary nonlinear field equations in terms of 
conjugate flow intrinsic coordinates is discussed in section [3] Section [4] deals with the existence 
of a principle of stationary action in the context of conjugate flow variations. Formal symmetry 
conditions and corresponding determining equations for the conjugate flow are finally obtained 
and discussed in section [5] for second-order nonlinear scalar field equations. 
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Figure 1. Sketch of the transport phenomenon of a volume of particles labeled by <t m 
according to two different realizations of the conjugate flow S M (a v ; u J ) , namely, x^(a"\ u^) 
and (cr v ; u J 2 ) ■ Shown are the trajectories of two different particles labeled as o\ and CT2- In 
this sketch, the time coordinate is not transformed, i.e. we have x° = a = t. 



2. The conjugate flow 

Conjugate flow is an intuitive physical notion which is represented mathematically by a 
continuous point transformation of (n + l)-dimensional (n denotes the number of spatial 
dimensions) Euclidean or Riemannian space into itself. In order to set up this transformation, let 
us consider a particle labeled by a v (y — 0, .., n, being the temporal component) and represent 
its trajectory in a fixed space-time Cartesian system as 

= x 11 {<j v ;u°) , /i,i/ = 0,..,n, (1) 

where u- 7 (j = 1, ...,N) is a vector field that solves a prescribed system of field equations. The 
transformation (JTJ) is assumed to be invertible (with diffcrentiable inverse) and eventually to 
possess even continuous derivatives up to a prescribed order except possibly at certain singular 
surfaces, curves or points. These requirements make (JTJ) a diffeomorphism, i.e. a time-dependent 
flow of curvilinear coordinates [351 HH 133 GUI 121] whose motion in space resembles in toto a 
physical fluid flow of classical mechanics. In figure [JJ we sketch two realizations of the conjugate 
flow (UJ for two different solutions fields corresponding, e.g., to different boundary or initial 
conditions. 

Coordinate flows being functionals of the solution to a field equation are obviously not 
new in the scientific literature. For instance, in the context of symmetry analysis of partial 
differential equations the so-called non-classical symmetries [5] are remarkable examples of 
solution-dependent transformations. Similarly, in classical Lagrangian fluid dynamics the 
trajectories of the fluid elements in space are obtained as local functionals of the velocity 
field IJi (x k ,t) that solves, e.g., the Navier-Stokes equations (x k is a fixed Cartesian coordinate 
system) 

dW V k^i^_^_ 1 d 2 W dU k _ () ._ i 

dt dx k dx k Re dx k dx kl dx k ' ' ' ' 

Such functional relation is defined by the solution to the well-known problem 



dX3{a\t) 
dt 



W (*V. *).*)> Xi(a\t ) = a j . (3) 
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Note that the physical fluid flow X 3 can be considered as a very particular type of conjugate 
flow since it is functionally dependent on the solution to Eq. @ by means of Eq. ([3]) . Another 
conjugate flow which is different from the physical fluid flow may be defined, e.g., by solving 

dX3 Q t l,t) = U j (x k (o-\t) + B"(o-\t),t) , X 3 (a\t )=a 3 . (4) 

where B u (a 1 , t) is a realization of a space-time Brownian motion. A remarkable result by Gomes 
[T7] shows that ensemble averaging - over random B" - of diffeomorphisms of type (jj) allows to 
construct a variational principle for the Navier-Stokes equations. 

2.1. Infinitesimal flow perturbations 

The components of the vector field u 3 appearing in Eq. (fT]) are, by definition, Cartesian 
components expressed in terms of conjugate flow intrinsic coordinates a v . In other words, if we 
denote by U 3 (x^ 1 ) the Cartesian components of the solution to a prescribed set of field equations 
expressed in Cartesian coordinates (e.g. Eqs. ©), then u 3 are defined as 

u j (a v ) d =U j (x"(a v )). (5) 

We remark that these are not the tensorial components [3111] of the vector field. Let us consider 
an infinitesimal perturbation of v? in the form 

u j (a v )+eip j {a u ), for e -> 0. (6) 

Then, to the first-order in e, we obtain the following perturbation in the conjugate flow 
trajectories (TTJ 

& {o v ; u 3 + etp 3 ) ~ 8" (<t v ; u 3 ) + e— tp> , (7) 

where, by definition 

S5P , def r x»(a»;u 3 +e(p 3 )-x»(a v ;u 3 ) , Q , 
— — w 3 = hm . (8) 

The quantity Sx^/Su 3 is known as Gateaux derivative [45 of the functional with respect to 
u? and under rather weak requirements (31) it is a continuous linear operator. The perturbed 
flow x^{a v \u 3 + eip 3 ) is schematically depicted in figure [2] and it is assumed to have the same 
regularity properties of the unperturbed one, i.e. invertibility and continuous derivatives up to 
prescribed order in all variables. 

Let us now postulate that the solution field u J is also functionally connected to the conjugate 
flow x^ and let us denote this functional relation by vP (a v ; x' 1 ). This fundamental assumption 
implies that an infinitesimal flow perturbation x^ + ecj)^ induces the following variation in the 
solution field u 3 

~ ■ Sn 3 ~ 

u J (tr";x^ + e^) ~ u 3 (a u ;x 3 ) + e— 0", (9) 

where, in analogy with Eq. ([8]). we have defined the Gateaux differential as 

5u3 r - f hm uJ{<7 "'^ + - ui ^ j ) (10) 

5x^ e->o e 

In the context of the Navier-Stokes equations, this means that a perturbation in the conjugate 
flow x 3 determines - by assumption - a perturbation in the velocity field U 3 that solves Eq. @. 
This ultimately results in a perturbation of the physical fluid flow X 3 by means of Eq. In 
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Figure 2. Conjugate flow perturbation + e0 M induced by a field perturbation + e<pi 
and corresponding deformation of the volume of particles u v transported by the flow. 



other words, by perturbing the conjugate flow in this case we are actually perturbing the physical 
fluid flow. At this point it is convenient to set 

? = S>< (id 

- St 11 

r = ^ (12) 

and write Eq. © and Eq. © as 

u j (a tl ;x v + e^ v ) ~ u j (a"; x v ) + e£» (cr"; 2") , (13) 
(a"; M 3 + e^') ~ x" (<t"; w J ) + (a v ; u 3 ) . (14) 

Note that in these equations we have </> M 7^ 0^ and (p J ^ (p ■ In fact, if we arbitrarily perform 
a simultaneous perturbation of u J and 3? M we cannot obviously expect that, in general, the 
functional disturbances arising from the Gateaux differentials (TTTT) and (|12|) coincide with the 
perturbations at the left hand side of Eqs. (Tl"3")) and (fT4"|). This immediately leads us to the 
question of which variable between v? and x 11 should be chosen as independent when performing 
perturbations. In the sequel we will be mostly concerned with perturbations induced in the 
conjugate flow through a variation of the solution field u J , i.e. we will mostly employ Eq. 
([Ml) , although the other approach, i.e. the one based Eq. (fT3|) . can be equivalently considered. 



3. Conjugate flow representation of field equations 



Several field equations of mathematical physics, remarkably the fluid mechanics equations, 
include naturally the concept of a conjugate flow within their formulation. Such a flow usually 
has a direct physical interpretation, e.g., trajectories of fluid elements in space, and it often 
constitutes the ground work on which dynamical results are constructed [I] [5]. Many other 
equations, however, do not include explicitly any term having a direct reference to a conjugate 
flow. This is the case, for example, of the classical heat equation, the Maxwell's equations of 
electrodynamics, the laws of elasticity and, undoubtedly, many others. The fundamental question 
at this point is whether it is possible to formulate a law that include both the field equations 
and the conjugate flow and allows to study their interaction, e.g., in the context of the principle 
of stationary action. The answer is affirmative and the simplest way to achieve this result is to 
represent the field equations relatively to a coordinate system which is advected by the conjugate 
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flow, namely, coordinates a v . In other words, we represent the partial differential equations 
governing the physical system on a curvilinear manifold which is assumed to be functionally 
dependent on their solution. 

This procedure obviously introduces a functional dependence in the form of the field 
equations from their solution. As a consequence, the equations look completely different in 
conjugate flow intrinsic coordinates and, in general, they are highly nonlinear. For example, by 
using the mathematical tools of | Appendix A[ it can be show that the classical one-dimensional 
heat equation 

dU d2U n (in 

-m- a ^ = > (15) 

where U(x,t) denotes the temperature in fixed Cartesian coordinates, can be written in terms 
of conjugate flow intrinsic coordinates as 

du 1 du dx a ( dx d 2 u d 2 x du\ 



dt dx/dadadt (dx/dof \da da 2 do 2 do) °' ^ 

where the flow a; is a functional of u. For illustration purposes, here we have assumed that 
the time variable is not transformed, i.e. we have set x° = a = t. By examining Eq. (|16[) 
under the conjugate flow perspective we see that a perturbation in the field u (a, t) induces also a 
perturbation in conjugate flow x(a, t;u) (see Eq. (fl4)l) and therefore the perturbed equation in 
conjugate flow intrinsic coordinates includes many terms arising from the perturbations of both 
u and x. Clearly, if the conjugate flow is in rest with respect to the fixed Cartesian coordinate 
system then Eq. (fTB")) coincides with Eq. (fT5|) . although the effects of the aforementioned 
functional perturbations are still present. 

The conjugate flow representation of a field equation is obviously much more complex than 
a standard formulation in fixed Cartesian coordinates. This has been observed, e.g., by Temam 
[36] . in the context of the Lagrangian representation of the Navier-Stokes equations. He pointed 
out that "the Lagrangian representation is not used too often because the Navier-Stokes equations 
in Lagrangian coordinates are highly nonlinear". Indeed, by using the results of | Appendix A| it 
can be shown that these equations can be written in general conjugate flow intrinsic coordinates 
as 

da- Ao+U da- Ak ~ da- + Re \ da^daP AkAk + da> daP Ak J ' {U > 

where the quantities A%, defined in Eq. (|A.6[) . are rather complex functionals of x^. Clearly, 
when the coordinate system a^ is advected exactly by the physical fluid flow, i.e. when the 
functional link between u 1 and x^ is defined by Eq. (J3J) , then Eq. (fT7|) coincides with the 
Lagrangian representation of the Navier-Stokes equations. 



3.1. Functional setting 

Let us associate with the physical system the linear function space U, whose elements are the 
iV-tuples u = (u 1 , .., u^) . Similarly, let us also consider the configuration space X, whose 
elements, denoted as x = represent (n + l)-dimensional conjugate flows, n being the 

number of spatial dimensions. In general, the configuration space is not a linear space because 
the summation of two conjugate flows is not a conjugate flow. This is due to the fact the 
superimposition of two invertiblc flows may not be invertible (the summation of two invertible 
Jacobian matrices is not necessarily invertible). However, the requirement that the perturbed 
conjugate flow has the same properties of the unperturbed one, i.e. it is still a diffcomorphism, is 
equivalent to state that locally, i.e. in the neighborhood of a particular flow x, the configuration 
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space X is linear or can be linearzed. In this sense we can say that the configuration space X is 
locally linear. Given this, an arbitrary field equation (or a system of field equations) written in 
terms of conjugate flow intrinsic coordinates can be synthesized as 

N £ (u) = v , (18) 

where Nj is, in general, a nonlinear operator while 0y denotes the null element of a third 
topological linear space V. The subscript x in reminds us that the operator is defined in 
terms of conjugate flow intrinsic coordinates a v , i.e. on the manifold x. For example, Eq. ((T6|) 
can be written in the form (fT5j) by defying Ng as 

def du 1 du dx a / dx d 2 u d 2 x du\ 

* (M) = dt ~ dxjfoda~~dt ~ (dx/da) 3 \da~da 2 ~ da^d^J ' ( ' 

The domain of the nonlinear operator Nj, is a space of functions satisfying the initial or the 
boundary conditions of the problem. In the conjugate flow theory, however, the operator N^(it) 
acts on both u and x and therefore it implicitly identifies two different domains, one within the 
space of fields hi and the other one within the configuration space X . These two domains will be 
denoted by D^(Ng) C hi and Dx{^x) C X, respectively (see figure[3]). The range of the operator 
Nj will be denoted by i?(Ng) C V. The next fundamental step in the functional setting of the 
conjugate flow theory of field equations is to introduce duality parings between the linear spaces 
hi, V and the locally linear one X through non-degenerate local bilinear forms [471 126) . To this 
end, let us define 

(v)„:Vx^I, (20) 
(■,■)£ :Vx X ->R. (21) 



The subscripts u and x in Eqs. (|20| and (|21j) emphasize the fact that such forms depend also on 
u and x, respectively (in a possibly nonlinear way). An explicit expression of (|20[) will be given 
in section l4~2l For the moment we simply observe that, locally, the forms (j2U)) and ([2~Tj) can be 
put in a correspondence through the linear transformations defined by Eqs. (|11[) and (|12[) . In 
fact, as shown in figure [31 the elements of -Dw(N^) in the neighborhood of a certain u are in 
correspondence with the elements of Dx(^x) in the neighborhood of a certain x. In practice, 
such a correspondence can be established locally through the linear operators Su/Sx and Sx/du. 
For instance, by using Eq. (|lip we obtain 

Su — — 

(v, (flu = (v, = (v, <t>)x- (22) 

We shall conclude this section by explaining why we have chosen the definition "conjugate flow" 
for the transformation (JTJ . To this end, we recall that the Gateaux differential of x with respect 
to u defines a linear functional from the space U to the space X = IA\ which is the conjugate 
space of U. Thus, for every admissible u G U, the flow x belongs to the conjugate space of 
U, hence the definition "conjugate flow". In a broader sense, the adjective "conjugate" simply 
emphasizes that there exist a functional relation between the flow x, the dynamic equations N of 
the field theory and their solution u. We also remark that a definition of conjugate flow already 
appeared in the scientific literature [3J [1] , as "a flow uniform in the direction of streaming which 
separately satisfy the hydrodynamical equations" . This definition, however, is clearly different 
from ours. 



3.2. Perturbation expansions 



By employing the operatorial approach developed in the previous section we can easily synthesize 
in a single operator equation the perturbative form of an arbitrary nonlinear field equation 
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Figure 3. Sketch of the function spaces employed for the functional setting of the conjugate 
flow theory. Shown are the domains D/j(N^) C hi and Dx{^s) Q X of the nonlinear operator 
representing the field equations. The range of Ng is denoted by R(N^) C V. We also show 
the correspondence between field perturbations (u + eip), conjugate flow perturbations (x + E(fi) 
and corresponding perturbations induced in the field equations (N_^ £ ^(m + £<p)) relatively to 
a specific representation (u, x, N^). The local bilinear forms that put the various spaces in 
duality are indicated in between the sets. 

(or a system of nonlinear equations) in the presence of a conjugate flow perturbation, i.e. a 
simultaneous perturbation of both the solution field and conjugate flow. To the first-order in e 
we have 

"<5I% <5N 5 



N 



X-\-£( 



(u + etp) = U S (u) + > 



Sx 



+ 



(23) 



where the Gateaux differentials appearing in Eq. (|2"3")) are defined as 

SN S def .. 1% (u + eip) - 1% (u) 
-ip = lim 



Su ' e->0 

-—0 = lim . 

ox e^O e 



(24) 



(25) 



provided that such limits exist. A function space representation of the conjugate flow 
perturbation is sketched in figure [31 



4. Conjugate flow action functionals 

Let us consider a field u G Du(N%) and a conjugate flow x e D;f(l%). The couple (u,x) does 
not necessarily have to be a solution to the field equation, i.e. Nj(u) ^ 0y- Disregarding the 
particular form of the operator 1%, it is useful to consider 

v = Ns(u) € i?(N £ ) (26) 

as a definition two vector fields [3T]: one in D^(Nj) and the other one D^(Nj), respectively. This 
allows us to introduce in a conceptually simple way the notion of a line integral of an operator 
according to a geometric standpoint which seems originally due to Volterra |48j . To this end, let 
us consider a one-parameter family of fields in the domain Du(N^) 

u = u x (0 < A < 1) . (27) 
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This can be regarded as a line in the function space Du(N^). With such line we can associate 
the number 

' (Hz(u x ),^) ax dA, (28) 

i.e. the path integral of the operator along the trajectory of functions u\ € Du(N%). We 
recall that (•, ■)„ in (|28[l denotes the local bilinear form ((20)) . In the context of the conjugate flow 
theory, we can also define the path integral of the operator along a trajectory of flows x\ in 
the space D ,*■(!%), i.e. 

h = J <Ns>(tt),^>2 x dA, (29) 

where (v)x denotes the local bilinear form (|2"Tj) . Therefore, we can define an infinitesimal 
circulation of an operator around a certain field u as well as around a certain conjugate flow x, 
these circulations being of course related by Eqs. (fTTj) and (fT"2"|) . If the line integrals (|2"B"|) and 
(|29p are independent of the path of integration then the operator N j is said to be potential with 
respect to the chosen local bilinear form. In this case the line integral from a prefixed element 
Mo to any element u in the domain of along an arbitrarily chosen path defines the action 
functional 



1 du^ 



A u [u] = A u [uq] + / (Nj? (u x ) , —)u x dX- (30) 



Similarly, the line integral from a prefixed conjugate flow xq to another flow x along an arbitrarily 
chosen line x\ defines another (dual) action functional 

A B [x] = A s [x ] + J\n Sx (u) , ^h,d\. (31) 

In turn, the operator Ng is said to be the gradient of the functionals v4 M [u] or This 
definition relies on the fact that if we calculate an infinitesimal line integral by using u or x as 
independent variable then we obtain, respectively, 

8A u [u] = {Nz(u),6u) u , SA s [x] = (N s (u),5x) s . (32) 

The relations f|32[) show that the equations of motion of the system, i.e. 1% (u) = 0y, can 
be obtained as as a stationary point of either or Ag[x], for arbitrary variations Su and 

5x, respectively. Thus, the theory of conjugate flows allow us to look for action functionals 
associated with field equations in two different ways, depending on which variable between u or 
x is assumed as independent. Clearly, if we consider u as independent then we are looking for the 
set of conjugate flows such that the field equation is potential. On the contrary, if we consider 
the conjugate flow x as independent then we are looking for the set of fields u such that the 
field equation is potential. In the sequel we will be mostly concerned with conjugate flows action 
functionals where the field u is considered as independent variable. 



4-1- Existence conditions 

In order to formulate the existence conditions of conjugate flow action functionals we follow the 
approach of Magri [2 6) . To this end, we consider two infinitesimal trajectories (two infinitesimal 
straight lines) of the field u in the function space Z?^(Ng) 

I : u — > u + tip 

II : u — > u + v^p 
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Correspondingly, we have the following infinitesimal conjugate flow perturbations 

I : x — > x + €(j) 

II : x — > x + vrj 

where <p and fj are related to the field perturbations ip and ip in the sense of Eq. (fl2l) . Due 
to this fundamental relation, an infinitesimal circulation of the operator around the element 
u G -Dw(Ns) is associated with an infinitesimal circulation of Ng around a flow x € Dx(N^). 
The vanishing of these simultaneous circulations with respect to the local bilinear form (PIUl) is 
synthesized by the condition 

(1% (u) , eip) u + (N £+e ^ (u + eip) , uip) a+elp = 

(1% (u) , vip) u + (N s+ „fj (u + vip) , etpju+vif,- (33) 
To the second-order in e and v we have 

Vj ^ # i im M^H^ (35) 



where 



denotes the Gateaux differential of the local bilinear form (j20j) . considered as a particular type 
of nonlinear operator on u. A substitution of Eq. (|3~if into Eq. (|3"3"f gives 

(^■<p + ip) u + (<p; Njf(u), V)« = 

oii da; 

(^■0+ ^rw)a + Nj(u),y>) u . (36) 
dii ox 

Finally, by using Eq. (JT2J) we can write the vanishing condition of the infinitesimal circulation 
entirely in terms of field perturbations tp and tp as 

(G s ip, ip) u + (ip; N s (u),ip) u = (G$ip, ip) u + (ip; N s (u),ip) u , (37) 

where the linear operator Gg is defined as 

G _ dcf SNs + SN$ $x / 38 % 
Su Sx 8u 

Thus, if the circulation vanishes along any infinitesimal closed line in Du(N%) then Eq. (|3"T)l 
must hold for every ip, tp and for all admissible u. This is the necessary condition for operators 
to be potential with respect to the local bilinear form f|20|) . If the domain of the operator is 
simply connected, then this condition is also sufficient. The happens, for example, when Djj(^x) 
is defined by linear homogeneous initial or boundary conditions (in this case Z^(IMg) is a convex 
set). We notice that Eq. (|3"T)) includes also interesting subcases. For example, if the flow x is 
not a functional of u then we have 

(nr^¥>! + (v; N 2i i>)u = (^"ip, <p)u + (tp; Ns, <p)u- (39) 

ou ou 
In addition, if the bilinear form (•, •)„ does not depend on u, i.e. if we are dealing with a global 
bilinear form, then the potential theory coincides with the classical theory of Vainberg [?S] . In 
this case the condition (|3"T)l reduces to 

(-7— ¥>,W = (-7— V,<P), ( 40 
ou ou 

namely, the Gateaux derivative of the operator N j must be symmetric with respect to the bilinear 
form (,). 
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4.2. The local bilinear form defining the conjugate flow action functional 

In classical, relativistic and quantum field theories the action functional has the standard form 
[221 1231 132] 

A = J C^da. (41) 

where C denotes a Lagrangian density, g is the determinant of the metric tensor and ^fgda is 
the invariant space-time volume element (da being a shorthand notation for da°da x ■ ■ -da n ). 
We recall that the square root of the metric tensor determinant is equal to the Jacobian 
determinant J of the transformation from fixed Cartesian to conjugate flow intrinsic coordinates 
(see |Appcndix A| for further details) . A comparison between Eq. (|4Tj) and Eq. (l30|) suggests 
that the local bilinear form to be considered for the conjugate flow formulation of the inverse 
problem of the calculus of variations is 

(a, b) u = f / abJda, aeU, beV (42) 



where the Jacobian determinant J is a rather complex functional of x (see Eq. (1A.5I0 . The 
domain £ appearing in the integral (|42[) is a four-dimensional volume of particles a v advected 
by the conjugate flow (see figure [1] and figure [5]). The form (|4"2"|) generalizes the bilinear form 
appearing in the Herivel-Lin variational principle [THl El UH EE) > where the volume of particles is 
advected precisely by the physical fluid flow. Note also that (l42j) is symmetric, non-degenerate 
and non- negative, i.e., it satisfies all the properties of an inner product. By using Eq. (|A.21[) we 
obtain the following Gateaux derivative 

def 

~de 



(<p;a,b) u = —[(a,b) u+eip } e=0 



abJV ■ 4>da (43) 



i2 



where, according to Eq. (|12[) . is a linear functional of (p, i.e. (|43l) is a trilinear form in a, b and 
ip. 



4-- 3. Incompressible flow perturbations 

Let us assume that the divergence of the perturbation field cf> appearing in Eq. (j4"3"]l vanishes. 
Under this assumption the symmetry condition (|37l) simplifies to 

(G s tp,ip) u = (G s ip,tp) u . (44) 

Equation (|44D basically requires the symmetry of relatively to the local inner product (14"2"1) . 
Thus, the application of the conjugate flow theory to the inverse problem of the calculus of 
variations is now reduced to look for an incompressible four-dimensional flow of curvilinear 
coordinates that symmetrizes the operator G^. 



5. Symmetrizing flows 

A field equation is said to be formally symmetric when the operator symmetry condition, e.g. 
Eq. (I44p . is satisfied disregarding the particular form of the boundary or the initial conditions 
associated with the problem. Clearly, when the domain of the operator is formed by a set of 
functions satisfying local homogeneous boundary and initial conditions then formal symmetry is a 
necessary condition for symmetry. Such a condition, however, is not sufficient even in the case of 
homogeneous boundaries |13) . In any case, it is useful to establish formal symmetry conditions 
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for particular classes of field equations. This has been done by Tonti [3HJ HO] using classical 
inner products in fixed Cartesian coordinates. In this section we obtain similar conditions in 
the context of conjugate flow variations. To this end, let us consider the following second-order 
nonlinear scalar field equation 

Ni(«) = /(M;w M ;v;^;^A) =°: ( 45 ) 

where the comma denotes partial differentiation with respect to a p (/i = 0,..,3), i.e. u ijL1 = 
du/dcr^. The Gateaux differential of Eq. (|45D is obtained as 

5N£ <5N£ 7 df df df df - u df , 

ar* + -*r* = ^ + du~/^ + duzr^ + wf v + a^f vX - {4<)) 



The conjugate flow perturbation <^ is related to the field perturbation <p by Eq. (TT2jl . We remark 
that, in general, such transformation could involve both derivatives and integrals. For example, 
it could be in the form 

' K p {o v ;u)^o + A" (u; u v ) p + («; a v ) <p, x . (47) 

The choice of the functional dependence between ^ and tp is actually a matter of investigation 
when looking for a conjugate flow variational principles. In this section we shall limit ourselves 
to algebraic flows, i.e. flows that can be expressed in the form 

x"=x*(o";u), (48) 

where the functions (to be determined) are local functionals of u, i.e. they do not involve 
integrals of the field u. Under these assumptions, the flow perturbation (f> p is easily obtained as 

At this point we set 



dr p 

" = ^-<p. (49) 
du 



def dx p 

du 



(50) 



def,,, . dK 

du 

This allows us to write the conjugate flow perturbation and its partial derivatives as 

= af*(p, (53) 
% = b£<p + a'><p, v , (54) 
fe, = < P V + Kv, P + W p <P,v + a^<p,u P . (55) 

Remarkably, the derivative order of cj>^ is the same as that of tp. In other words, the fcth-order 
derivative of u with respect to a v involves the fcth-order derivative of x^ (see Eqs. (|A.14|) and 
(|A.15|I ). This is why we have included the second-order derivative of the conjugate flow in the 
second-order scalar field equation (|4"5)) . The symmetry condition (|4"4"|) can be now explicitly 
written as 

/ ipHtpJda + / ipB^tp^Jda + / ipF plJ tp. pv Jda = 
Js is Js 

/ pHiPJda+ { pB p ifj :P _Jda + j tpF pv ^^ pu Jda, (56) 
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where H, £? M and F pv are obtained as 

H - du + dx^ + dx%°^ (57) 

Next, we integrate by parts the integrals at right hand side of Eq. (|56"]l and we neglect all the 
boundary terms (we are looking for formal symmetry). By using the identity (fP2|t we obtain, 
for instance 

I (pB^^Jda = - I tpB u ip^Jdcr^- I ipB^ipJda- j i/jBT^Jda, (60) 

where denotes the affine connection of the conjugate flow. Proceeding similarly with the 
other terms at the right hand of Eq. (156|) . we conclude that the formal symmetry requirement 
is satisfied if and only if 

F vp = F p \ (61) 
B" = jW + r^F"' (62) 

s£ = + 2r^Fj + r^r^F^ + r^ PtV F v " - t^b». (63) 

A differentiation of Eq. (162)) with respect to er M and subsequent substitution in Eq. (IBlfl) yields 
the single relation 

(F%>+rl p F»»-B»)rZ I/ = 0, (64) 
which is equivalent to the following system of determining equations for the conjugate flow 

Kp +V Pp FPV = B " ■ ( 65 ) 
The system (|65p defines the relation x IM (cr l ';u), i.e. the Lie group manifold that guarantees the 
existence of a principle of stationary action for the field equation (|4"5|) . Once x^icr 1 '; u) is available, 
the action functional of the field theory can be explicitly constructed by calculating the integral 
(f3T)]) along an arbitrary trajectory of admissible functions u\. This yields an action of type (|4"Tj). 

We notice that if we remove the functional link between x^ and u, then the conditions ([55]) 
consistently reduce to those of Tonti |39[ |4"0] . In order to see this, we simply set a M , b% equal 
to zero in Eq. (|35|) and Eq. {55]) and then substitute them into Eq. ([6"5]) . The result in fixed 
Cartesian coordinates (r^ p = 0) is 

J_ ( df \ _df_ =Q (66) 
dx p \du^ v J du, v 

This is the classical condition arising from the symmetry requirement of a second-order nonlinear 
scalar equation [131 US] . 



6. Summary 

We have developed a new approach to construct an action functional for a field theory described 
in terms of nonlinear partial differential equations (PDEs). The key idea relies on an intrinsic 
representation of the PDEs governing the physical system relatively to a diffeomorphic flow of 
coordinates (the conjugate flow) which is assumed to be a functional of their solution. This flow 
can be selected in order to symmetrize the Gateaux derivative of the field equations relatively to 
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a suitable (advective) bilinear form. This is equivalent to require that the equations of motion 
of the field theory can be derived from a principle of stationary action on a Lie group manifold. 
By using a general operator framework, we have obtained the determining equations of such 
symmetrizing manifold for a second-order nonlinear scalar field theory and shown that they are 
consistent with classical results in fixed Cartesian coordinates. Once the symmetrizing manifold 
is available, then the conjugate flow action functional can be constructed explicitly through path 
integration. In particular, the duality principle between the conjugate flow and the solution field 
discussed in section U] allows us to perform integrations either in terms of flows or in terms of 
fields. 

The proposed new methodology can be generalized to vectorial and tensorial field theories. 
In particular, it can be applied to the Navier-Stokes equations, for which a great research effort 
has focused in obtaining a physically meaningful principle of stationary action [5] [UJ HS1 HI] ■ 
Recent results of Gomes [T71 [TS] , Eyink [12] and Constantin [TUl E] indeed have shown that an 
action principle can be constructed for the Navier-Stokes equations on random diffeomorphisms 
[TBI l50l [33] . These random flows are usually defined in terms of perturbations of a Lagrangian 
base flow. In the proposed new framework, the variational principle for the Navier-Stokes 
equations may be constructed on a generalized diffeomorphism (Lie group manifold) satisfying 
a set of PDEs similar to (R)5l) . These equations, however, are strongly nonlinear and their study 
will be the objective of a future work. 
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Appendix A. Representation of field equations in conjugate flow intrinsic 
coordinates 

We recall some fundamental identities from differential geometry [331 E US HI] that allow us to 
write the dynamic equations of a physical system in terms of conjugate flow intrinsic coordinates. 
To this end, let us first consider the Jacobian of the conjugate flow transformation (JTJ) 

(A.1) 

It is easy to verify that the transpose of the algebraic complement of Jjj has tensorial expression 
(repeated indices are summed) 

r x dof 1 w... d ^ 9x X ( a ^ 

C p - -e z pilX ..——..., (A.2) 

where n denotes the total number of spatial dimensions while e Xva " and e p ^\... are multi- 
dimensional permutation symbols, i.e. Levi-Civita tensorial densities. In particular, if we 
consider only two dimensions (e.g., one spatial and one temporal dimension) then we obtain 
the simple expression (all indices are from to 1) 

C p =e "ep^T^j- (A. 3) 

Similarly, in 1 + 3 dimensions, i.e. one temporal and three spatial dimensions (all indices are 
from to 3) 

n \ _ l \u a B dx x dx s 
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By using Eqs. (|A.1[) and (|A.2[) we obtain the Jacobian determinant 

jH f det(J^) = ^ I J^. 
This allows us to write the the inverse of the Jacobian matrix (IA.1I) as 



A x def Cp_ 
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(A.5) 



(A.6) 



o J 
We denote by 

a" =a" (x";u j ) , n,v = 0,..,n (A.7) 

the inverse transformation of ([1]). Such inverse transformation exists and is differentiable by 
definition of conjugate flow. From the well known identity 
&a v dr* 

s?£*-« .... (A : 8) 

we obtain the following fundamental expression of the partial derivatives da" /dx* as a function 
of a x 

(A.9) 



where A 1 ^ is defined in Eq. (|A.6j) . It is useful to write down A"^ explicitly for the two-dimensional 

dx x /do x -dx°/da x ] 

-dx x /da° dxP/da° \ ' (A " iUj 
dx 1 dx° dx° dx x 



A° A\ 
A 1 A 1 



where 



J 



da 1 da da 1 da ' 
If time is not transformed, i.e. if x° = a° = t, then Eq. (|A.10[) reduces to 

1 

-(dx/dt)/(dx/da) \/{dx/da) 

where we have denoted by a = a 1 and x = x 1 . 



\A° 


A\ - 




. Ah 


A\ _ 





(A.ll) 
(A.12) 



Partial differentiation in intrinsic coordinates 

We consider a vector field [P (x^) expressed relatively to a fixed Cartesian coordinate system x*. 
If we express the a^-dependence of U 3 in terms of the trajectories of the particles a" (advected 
by the conjugate flow) we obtain the following equivalent representations 

U j (x") = U j (£" {a")) = u ° = u3 (v U («'*)) • ( A -!3) 
The transformation law for partial derivatives of U 3 ' is obtained by differentiating Eq. (|A.13|) 

dx* da" dx* da" V ' ' 

where the quantities da" /dx* are expressed in coordinates a" through the fundamental relation 
(|A.9[) . Let us now evaluate the second derivative with respect to x v and express the result in 
conjugate flow intrinsic coordinates. To this end let us perform an additional differentiation of 
(IA.14I) with respect to x". This yields 

d 2 U j d 2 u j 

" x "- "'" (A.15) 



ap A x 

dx*dx" da x daP v " ' da x da* 

By using the expressions of J and C"^ obtained in (|A.5|) and (IA.2[) it is possible to manipulate 

Eq. (IA.15[) further. However, it is more convenient to obtain first A" explicitly as a function of 

a* and then perform the differentiation appearing in (|A.15[) . 
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Perturbations in the metric tensor, affine connection and Jacobian determinant 

When the conjugate flow ^ undergoes an infinitesimal disturbance of type (|14D then all the 
quantities related to its intrinsic geometry are subject to small variations. For instance, the 
metric tensor 

def dx? dx 13 

becomes, to the first-order in e, + eh pu where 



def 



The corresponding perturbation in the affine connection (Christoffcl symbols of the second kind) 

def 1 ap f dg pu dg pv _ dg uv 
2 9 V do v da^ do-P 



r£„ = =;9 ap ( S + ir£ ~ S ) (A. 18) 



is 

ST^ - -eg ?h p ^ v + -g ^— + — - — j . (A.19) 
This can be also expressed in a covariant form as [33] 

ST%, = \g ap (h ptlv , + V;m " W) ( A - 2 °) 

the covariant derivatives ";" being of course constructed by using the unperturbed affine 
connection T" . These results allow to compute the conjugate flow perturbation of other 
fundamental geometric quantities such as the Riemann- Christoffcl curvature tensor. Next, we 
determine the perturbation of the Jacobian determinant (| A.5|) induced by a small disturbance 
in the conjugate flow. To this end, we substitute Eq. (fT?| into Eq. (|A.5[) and we keep only the 
terms that are linear in e. This yields 




where def)^ jdx^ denotes the divergence of the flow perturbation (remember that ^ are Cartesian 
components). Another useful formula involving the Jacobian determinant is 

g£ = JI V (A.22) 
This can be easily proved by directly differentiating Eq. (|A.5[) with respect to <t m 

da^ p do v do» dxP da v da^ ' \ ■ ) 
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